Constraints on dark matter annihilation by radio observations of M31 
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We used radio observations of the neighbour galaxy M31 in order to put constraints on dark matter 
particle mass and annihilation cross section. Dark matter annihilation in M31 halo produces highly 
energetic leptons, which emit synchrotron radiation on radio frequencies in the galactic magnetic 
field. We predicted expected radio fluxes for the two annihilation channels: XX ~^ ^^ ^ind XX ~^ 
T^T~ . We then compared them with available data on the central radio emission of M31 as observed 
by four radio surveys: VLSS (74 MHz), WENSS (325 MHz), NVSS (1400 MHz) and GB6 (4850 
MHz). Assuming a standard NFW dark matter density profile and a conservative magnetic field 
distribution inside the Andromeda galaxy, we find that the thermal relic annihilation cross section or 
higher {av) > 3 • 10~^^ cm'^/s are only allowed for WIMP masses greater than ^ 100 GeV and ~ 55 
GeV for annihilation into bb and r~^r~ respectively. Taking into account potential uncertainties in 
the distributions of DM density and magnetic field, the mentioned WIMP limiting masses can be as 
low as 23 GeV for both channels, and as high as 280 and 130 GeV for annihilation into bb and r~^r~ 
respectively. These mass values exceed the best up-to-day known constraints from Fermi gamma 
observations: 40 GeV and 19 GeV respectively [A.Geringer-Sameth and S.M.Koushiappas, Phys. 
Rev. Lett. 107, 241303 (2011)]. Precise measurements of the magnetic field in the relevant region 
and better reconstruction of the DM density profile of M31 will be able to reduce the uncertainties 
of our exclusion limits. 

PACS numbers: 95.35.+d, 98.56.Ne, 95.85.Bh 



I. INTRODUCTION 

The physical nature of the dark matter (DM) is contin- 
uing to be not understood. The most probable candidate 
for the role of DM is a Weakly Interacting Massive Par- 
ticle (WIMP): a supersymmetric partner of a Standard 
Model (SM) particle. There are three main approaches in 
attempts to detect WIMPs: accelerator searches, direct 
searches and indirect searches. Although great efforts 
have been dedicated already to all of these directions, 
the DM puzzle still seems to be very far from its final 
and clear solution. 

This article is related to the last mentioned search 
strategy: indirect detection. The idea of indirect detec- 
tion is based on the opportunity of pair annihilation of 
WIMPs causing production of different highly energetic 
SM particles, which eventually decay into stable parti- 
cles like leptons, protons and others. These yields may 
emit electromagnetic radiation through various mecha- 
nisms like synchrotron emission. Inverse Compton Scat- 
tering (ICS), etc. in different astrophysical objects. By 
detecting such radiation we can infer WIMP properties. 
Since there are evidences of DM existence in any ob- 
ject from dwarf galaxies to largest galaxy clusters, any 
of these objects is potentially a good target for indi- 
rect DM searches. Indirect searches of dark matter have 
been already extensively exploited considering all wave- 
length ranges from radio to gamma rays (see [1] for a 



review). Conventional ACDM model of cosmology sug- 
gests WIMP velocity averaged annihilation cross section 
(av) ^ 3 • 10~^^ cm^/s, which yields the correct cur- 
rent DM abundance in the Universe in case the DM is 
thermally produced (see e.g. [1]). Indirect searches may 
either lead to the discovery of DM with given particle 
annihilation cross section and mass ^, or constrain the 
cross section to be below the level of a thermal relic over 
all plausible range of WIMP masses (from several GeV 
to several TeV). The latter case may imply that in fact 
we do not understand the DM phenomenon in cosmology 
all that well, and a more sophisticated DM production 
mechanism is needed. 

Currently the most promising direction in indirect de- 
tection seems to be related to Fermi searches of the pri- 
mary gamma emission produced directly by WIMP an- 
nihilation (see [3]). At the moment, these constraints are 
the strongest: they exclude WIMP masses smaller than 
about 40 GeV and 19 GeV for bb and r~^r~ annihila- 
tion channels respectively for the standard thermal relic 
annihilation cross section. While the allowed region of 
WIMP parameters space is gradually shrinking, WIMP 
masses of hundreds and thousands GeV are still allowed. 
Indirect searches may allow further insights in the DM 
phenomenon. Wavelength bands ranging from the radio 
to the gamma rays and various astrophysical objects may 
be exploited further to this aim. In this work we studied. 
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in particular, radio observations of the closest big galaxy 
M31. Relativistic leptons produced by DM annihilation 
in its halo emit synchrotron radiation on radio frequen- 
cies due to Andromeda's magnetic field. Given the size 
of the galaxy and its DM halo, the strength of its mag- 
netic field and its close distance, large radio fluxes from 
DM annihilation may be expected from M31. There- 
fore MSl's radio properties may be exploited in order 
to put strong upper limits on the annihilation cross sec- 
tion. Surprisingly, almost no attempts in this direction 
have been made in the past: in the literature only few 
articles can be found dedicated to indirect searches in 
M31. The most remarkable among them is [4]. These 
authors obtained some constraints on WIMP parameter 
space by Cherenkov ground based gamma observations. 
However, those observations did not have enough sen- 
sitivity to probe the relevant region of DM parameters, 
and since then no significant progress has been made on 
this object. 

We computed the expected radio flux due to DM anni- 
hilation (section II), and then compared it with available 
data of radio observations of M31, which allowed us to 
put upper limits on the annihilation cross section (section 
III). For comparison with real observations we chose all 
appropriate radio surveys, which cover a wide range of 
frequencies: VLSS (74 MHz), WENSS (325 MHz), NVSS 
(1400 MHz) and GB6 (4850 MHz). The limits obtained 
can be considered as conservative because we did not 
make any specific assumptions about the radio emission 
other than the one from DM in the center of M31, and 
allowed for an unconstrained contribution from all other 
unknown backgrounds. 

We calculated the constraints for two annihilation 
channels, specifically the ones annihilating into bb and 
r~^r~ pairs. We chose these channels among all possibil- 
ities because, as explained e.g. in [5], bb and r~^r~ nearly 
present the channels with the softest and hardest lepton 
yields respectively. Any other case would therefore pro- 
duce radio fluxes at intermediate levels with respect to 
these two. In this sense, these can be considered as the 
two limiting cases. 

In our analysis we considered only the central part of 
M31, and specifically the bulge area of circular shape 
with the angular radius a ^ 5' around galactic center (see 
fig. 1). We chose this specific region of interest (ROI) 
as a target of indirect searches because of the following 
considerations: i) the radio quietness of the M31 nucleus, 
which indicates low contamination in the radio band by 
other standard astrophysical processes; ii) the absence of 
any projected point source inside it; in) sufficient halo 
size to produce a relevant signal. More details about the 
ROI choice will be explained below. 

Through our paper we adopted the Hubble constant 
value Hq = 71 km/(s-Mpc), which was taken from the 
WMAP7 data at http://lambda.gsfc.nasa.gov/. 

This paper is organized as follows: section II describes 



computation of the expected radio emission properties 
due to WIMP annihilation in M31, in section III we de- 




Figure 1. Optical image of M31 with the marked region of 
interest (ROI), selected for our purposes. Circle radius is 5^ 



rived the actual constraints by comparison of the pre- 
dicted fluxes with observational data, and section IV 
summarizes the results of our work. 



II. COMPUTING THE RADIO FLUX 
A. General theory 

In this section we presented the procedure for the com- 
putation of the radio flux density from the center of M31. 
We neglected here the potential absorption of radio emis- 
sion between the source and the observer, since our es- 
timates showed that it occurs at a negligible level (see 
appendix A). In the case of an optically transparent emit- 
ting medium, the total flux density from our ROI in M31 
can be obtained just by integrating the local medium 
emissivity j(z^, r) over the volume of halo contained in 
our ROI: 



S = 



j j{v,r)dV 



(1) 



where d = 785 ±25 kpc is the distance between us and the 
Andromeda center ([6]), r is the position vector inside the 
M31 halo originating in the M31 center. We disregarded 
here all (small) redshift effects. Then as a next step we 
needed to compute the local emission coefficient at an 
arbitrary position in M31 halo jiy^f). The synchrotron 
emissivity of leptons produced by WIMP annihilation has 
the form: 
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where e^ and e~ represent positrons and electrons re- 
spectively, Pe{E^v^r) is the synchrotron emission power 
of one lepton with energy £^ on a frequency v (measured 
in erg/(s-Hz) in CGS), ^^ is the energy distribution of 
leptons - the number of leptons per unit volume per unit 
energy range. We assumed electron and positron terms 
in eq. (2) to be equal to each other, which is reasonable 
because these both species behave similarly in all relevant 
aspects. According to [7]: 
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;|^ — plasma frequency of ambient plasma, 

^ 1.25(648)-^/-^^t-^/^e~^ - synchrotron spectrum, 

-^ — lepton Lorentz-factor. 

(4) 
This set of formulas describes the synchrotron power 
of a lepton in magnetic field P(r) in the presence of an 
ambient plasma with concentration n. Integration over 
0' in eq. (3) represents the averaging over all possible 
random angles between lepton's velocity and magnetic 
field. In order to proceed our computation of the radio 
flux expected we needed to derive the energy distribution 
^^. For this purpose we applied the standard diffusion 
equation: 

d dUe ^ ( ^^ dUe \ d f ^ ,^ . dUe \ / ^ X 

(5) 
where D is the spatial diffusion coefficient, b{E^ f) is the 
energy loss rate for a lepton through various energy dissi- 
pation mechanisms (measured in the units of energy per 
time), qe{E^r) is the source function - how many elec- 
trons (or positrons - but only one of these two species) 
are produced by WIMP annihilation per unit time per 
unit volume per unit energy range. We can simplify this 
equation assuming the stationary limit ^ =0, as com- 
monly done for diffusion of DM annihilation products in 
galaxies and clusters (see e.g. [8]). Another useful simpli- 
fication is the absence of the spatial diffusion of annihila- 
tion products, which will make the first term in the r.h.s. 
vanish. We investigated the validity of such assumption 
by comparing the characteristic diffusion length of newly 



injected leptons with the size of our ROI and concluded 
that it could be done without significant effect on our 
results (see appendix B for details). With these two sim- 
plifications we can easily solve the diffusion equation: 
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(6) 



At this point we should specify the functions b{E^f) 
and qe{E^r). The energy loss rate b{E,f) is constituted 
mainly by four different cooling processes: ICS emission, 
synchrotron emission, bremsstrahlung and Coulomb scat- 
tering: 



b{E,r) = bics^bs 



Obre 



-bcoi 



(7) 



Using previous results presented in the literature ([7], [1] 
and [9]) we constructed each of these four terms respec- 
tively as follows: 



bics = 0.76 



Ur, 



1 eV-cm 



1 GeVj ' 



(8) 



here and in other loss terms numerical pre-factors like 
0.76 follows from the exact analytical computation of the 
corresponding loss rate for a lepton. Uph denotes the 
total energy density of radiation at relevant locations. 
For the galactic center region it is constituted mainly by 
star light photons. We substituted for this quantity the 
fixed value Uph = 8 eV-cm~^, which is quoted in [1] as a 
characteristic value for the MW center. 
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(10) 



where n denotes the concentration of ambient plasma 
in cm~^.^ For this quantity we also used the constant 
value n = 0.1 cm~^ relying on the results of study [10], 
where such plasma concentration was derived as a typical 
average value for the Galactic center region. 
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^ This loss term 6{,^em was presented in another form in the work 
[7]. However, we suspect a mistake at this point in [7], because 
our verification of this loss rate did not confirm the expression 
there. For this reason, we use here the expression precisely de- 
rived in [9] (formula (6.74) there), which appears to have better 
justification. 



All loss rates here are measured in the units 10~^^ GeV/s. 
For illustrative purposes, we showed in fig. 2 the energy 
dependences of each term in eq. (8)-(ll) over the rele- 
vant range of lepton energies. We can clearly see that the 
synchrotron losses term depends on the location through 
B{r). For this reason, we presented bgync on fig. 2 for 
three relevant trial locations, which are away from the 
M31 center by 0, 1 and 2 kpc. More details about B{f) 
distribution will be discussed in subsection II C. In the 
most general case, the quantities Uph and n depend on 
f as well, indeed. However, at the current level of accu- 
racy of our model we made two simplifying assumptions. 
First, we assumed constant values for these quantities 
with conservative choices; and second, we adopted the 
values for the MW relying on high similarity between 
these two galaxies. 

As for the source function qe{E^r)^ it is computed es- 
sentially as a product of number of leptons produced by 
one WIMP annihilation per unit energy range ^^ and 
number of annihilations per unit time per unit volume: 



qe{E,r) = ^(^ 



pDM(r) 



m^ 



■ .dNe 



(12) 



where Pdm{^) is the local DM density, {av) is thermally 
averaged annihilation cross section and m^ is WIMP 
mass. Annihilation yields of leptons ^^ for two neces- 
sary channels XX ~^ ^^i t^t~ were taken from the web- 
site http : //www . marcocirelli . net/PPPC4DMID . html. 
A detailed description of this resource is presented in 
[11]. We took the pre-computed annihilation yields from 
this resource in their last version, which includes elec- 
troweak corrections. The importance of these corrections 
for a yields computation was justified in [12]. Inclusion 
of electroweak corrections leads to slightly different sec- 
ondary lepton yields from annihilation into primary prod- 
ucts with different polarizations: left- or right-polarized 
tau- leptons, transversely- or longitudinally-polarized W- 
bosons, etc. But as practice showed, the difference in the 
final radiation fiux is very minor due to this splitting. 
For example, the final fiux for the case of M31 DM halo 
concentration cioo = 12, observational frequency v = 7A 
MHz, WIMP mass m^ = 100 GeV and annihilation into 
T^T~ differs just by < 1% for different polarizations; 
which is much smaller than other total fiux uncertainties 
presented. That's why further we did not distinguish 
different polarizations of primary annihilation products. 
For the next step we needed to specify the DM den- 
sity distribution Pdm{t)^ which appears in formula (12). 
Next subsection is dedicated to the DM density profile 
derivation. 



B. DM density distribution 

As a model of DM density distribution Pdm{t) in M31 
halo we used the standard NFW profile, firstly introduced 
in [13] and widely used since that due to its universality. 
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Figure 2. The energy dependence of various types of lepton 
energy losses. Synchrotron losses are shown for three repre- 
sentative locations around the M31 center. For more details 
see subsections II A and II C. Color version of our article is 
available online. 



Some other profiles potentially describe a DM distribu- 
tion in galactic halos even better than NFW - this was 
demonstrated e.g. in [14] for the Einasto density pro- 
file. However, we did not use any other profiles except 
NFW further due to the following reasons. First, all cur- 
rently available models of the DM density distribution 
in M31 were made in the frame of NFW profile. The 
second reason is that among all profiles, the NFW yields 
moderate emission fiuxes due to DM annihilation (lower 
than e.g. Moore profile - see [8]), and this makes our 
constraints conservative with respect to the assumptions 
on the adopted density profile. Also, the NFW profile is 
specified by only two free parameters, while a larger num- 
ber of parameters is required by some other profiles. This 
choice therefore minimizes the overall number of "degrees 
of freedom" in our model. 

Quantitatively the DM density distribution of our 
choice is the following: 



Ar) 



Pdm{t = 50 pc), if ^ r < 50 pc; 
^ if r ^ 50 pc. 






(13) 



It presents the exact NFW profile with a minor modifica- 
tion, which can be called as a fiat core: we do not extrap- 
olate the density below the radial distance of rt = 50 pc, 
inside this central region we leave DM density at the con- 
stant level of pDAdi^ = 50 pc). This choice is motivated 
by the divergence of the NFW distribution towards the 
central point of the halo, which prevents to estimate the 
DM density reliably in the very central region. The NFW 



profile provides reliable density estimate down to the ra- 
dius ~ 50 pc for a M31-size halo (see e.g. fig. 4 in [15]). 
That's why in order to be conservative in our computa- 
tion of the medium emissivity, which depends quadrati- 
cally on the local DM density, we truncated NFW profile 
at r = rt = 50 pc and left the density at smaller radii 
on the level of Pdm{t = 50 pc). We tested the depen- 
dence of our results on the specific arbitrary choice of the 
truncation radius Vf. This study showed that a decrease 
(increase) of rt by a factor two (a factor four) produce 
an increase (decrease) of ^ 3% (~ 10%) in final fluxes, 
which is far below the overall level of accuracy of our 
model. 

The DM density profile eq. (13) is completely defined 
for a specific halo by two parameters - scaling radius r^ 
and scaling density pg. Thus, we needed to estimate 
them for M31. These two parameters are unambigu- 
ously linked to another two halo parameters, which are 
more meaningful - the halo mass Ma and its concentra- 
tion ca- a here denotes an overdensity: halo mass Ma 
by definition means the DM halo mass enclosed in the 
sphere, average density inside which is equal to Apcrit 

with pcrit = -g:;^ = 9-5 • 10~^^ g-cm~^ being the current 
critical density of the Universe. In our article we used 
A = 100 due to some practical circumstances. How- 
ever, the parameters of the actual DM density distribu- 
tion (13) do not depend indeed on a A choice. Halo mass 
and concentration are connected with the NFW profile 
parameters through the following relations: 
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Ca = — , Ma = / pDM{r)47rr dr, 

(14) 
As for the determination of the relevant parameters for 
M31, we used the results of Mioo and cioo from [6] and 
[16]. While these are in a good agreement with each 
other in determination of Mioo ((1.2 ± 0.3) • lO^^M© vs. 
(0.91 ± 0.16) • lO^^M© respectively), they show a rather 
big discrepancy in the estimation of cioo (the most prob- 
able values cited being ^ 12 and ~ 28 respectively). The 
expected radiation fluxes from these two concentration 
values differ by about one order of magnitude. This is 
why we decided to treat these two cases separately. They 
can be considered as two limiting cases yielding the most 
conservative constraints for cioo = 12 and the most opti- 
mistic constraints for cioo = 28. 

At this point we were ready to specify the NFW profile 
by estimation of its two parameters Ps , r^ . For the first 
case Cioo = 12 we took Mioo = 1.2 • lO^^M© from [6] and 
obtained Ps,rs by formulas (14). For the second case 
Cioo = 28 we did not need to calculate ps^Tg - they were 
taken directly from the table 2 in [16] as the best-fitting 
values. 

Another potentially relevant question in this subsec- 
tion is substructures contribution in our flux. It is well 
known that any DM halo contains a lot of small subhalos 
inside it (see e.g. [1]). These subhalos are very dense 



and numerous, that's why they are able to substantially 
increase the total flux due to DM annihilation from a 
whole halo - by 10-100 times (see e.g. [17]). But in our 
case we can neglect substructure contribution, because 
the main part of our expected flux comes from very cen- 
tral region with a size about 1 kpc. As can be seen e.g. 
in [17], substructures do not survive so close to the cen- 
ter due to tidal disruption, and their contribution to the 
total expected flux would be negligible at so small radii. 
Thus we did not need to include subhalos in our anal- 
ysis. Moreover, subhalos' presence would only increase 
the overall flux, so that the constraints obtained here are 
to be considered conservative. 

Thus, at this point we have completely specified the 
DM density distribution in M31 halo. Now we can move 
to the next step - specifying the magnetic field distribu- 
tion. 



C. Magnetic field 

The emission due to DM annihilation, which we hope 
to detect, is generated through the synchrotron mecha- 
nism. For synchrotron emission modeling it's crucially 
important to know the magnetic field strength distri- 
bution in our emitting volume. A global axis symme- 
try of the large spiral galaxy M31 naturally suggests as- 
sumption of two-parametric magnetic field distribution 
B{p'^z')^ where B depends on the radial distance from 
the center p' in the galactic plane and the vertical height 
z' above the galactic plane. 

As for B dependence on the vertical coordinate z' ^ 
it's plausible to assume exponential dependence B ~ 
exp(— |z'|/zq) with some scale height Zq. Such depen- 
dence is commonly used for large spiral galaxies - see 
e.g. [18]. In order to consider a reasonable radial de- 
pendence of B in the galactic plane we used the findings 
of [19], [20]. In [19] large-scale magnetic field in the disk 
plane was measured between radial distances of 6 and 14 
kpc. These measurements yielded almost constant mag- 
netic field over the whole mentioned annular region of 
^ 7 /iG, which is the expected value for galactic discs (see 
e.g. [21]). As for the magnetic field strength in the inner 
region p' < Q kpc, we were unable to find reliable data. 
Some measurements for this region were reported only in 
[22] and only for two trial locations with poor justifica- 
tion. Besides this information on the central field prop- 
erties we can reliably assume probably only one thing - 
that the field grows towards the center. In such situa- 
tion we made a decision to introduce the following field 
distribution (including the vertical dependence part): 



B{p',z')=[B 
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(15) 
Here Bio = 7 pG - the plateau value observed on 

p' ^ 6 — 14 kpc. The term ^^expf— -^j presents 

an exponentially decreasing central spike with the char- 



acteristic radial extent Pq. Of course, such distribu- 
tion has a non-physical plateau extending to infinity 
B{p' = oo^z^ = 0) = 7 fiG. But this plateau does not 
matter for our calculations: at distances of tens of kpc 
DM density becomes very small, and a contribution of 
these regions to the overall flux is tiny. The distribu- 
tion in eq. (15) can be completely specified by two 
free parameters: the vertical scale height Zq and the 
central field value 5(0,0) = Biq -\- Eg. After speci- 
fying these two parameters, the radial scale length p^ 
becomes automatically set by smooth connection of the 

central exponential spike Bg exp ( — A- ) to the plateau 

B{p' >6kpc,^' = 0) ^7/iG. 

Now let's discuss the specification of the values of Zq 
and 5(0,0). We estimated the scale height Zq using [23] 
and [19]. The first mentioned article outlines such a gen- 
eral property of galactic magnetic fields: the scale height 
Zq is typically about 4 times greater than the scale height 
of a synchrotron emission, which arises from cosmic rays 
in a whole galaxy volume. The last mentioned syn- 
chrotron scale height was reported in [19] for M31 with a 
value ~ 0.3 kpc. This would yield Zq ~ 0.3 kpc • 4 = 1.2 
kpc. Moreover, such a value for Zq is confirmed by the 
authors of [19] on the basis of other independent con- 
siderations: they conclude that the minimal expected 
Zq value is not less than ^ 1 kpc. We therefore chose 
z'q = 1.2 kpc for all our calculations. This estimate re- 
flects a conservative choice, since lower values seem to be 
unrealistic and possible greater values can only increase 
the final expected flux and, hence, strengthen the final 
constraints. 

As for the second necessary parameter - 5(0,0), the 
situation is much less certain. As we already men- 
tioned, we found in the literature only one attempt to 
measure the central field values in M31 - [22]. This 
work reports the total field strengths for the two lo- 
cations: B{p^ ^ 0.3 kpc, z^ ^ 0) = 15 ± 3 fiG and 
B{p' ^ 0.9 kpc, z' ^0) = 19±3 /nG. These values were 
obtained assuming energy equipartition between cosmic 
ray particles and magnetic fields. Such assumption, how- 
ever, has been shown to lead to an underestimation of the 
magnetic field in similar setting. For example, the au- 
thors of [20] explored the central field properties of the 
MW and derived the stringent lower limit of 50 /nG on 
the field strength based on multi- wavelength studies of a 
cosmic ray emission. At the same time, they mentioned 
that a simple equipartition assumption implies a central 
field strength of only ~ 10 /iG. As for the upper bound 
on a central field strength, [20] demonstrates that there 
is no certainly known value for it. Values of few hundreds 
/iG are mentioned to be absolutely possible. 

Considering the possible underestimation of M31 mag- 
netic field due to the equipartition assumption ([22]), 
and the similarity between MW and M31, we chose 
5(0, 0) = 50 /iG as the most probable value. However, 
we performed all our calculations for alternative field val- 
ues as well (which are considered as less probable) in or- 
der to study the dependence of our final results on this 
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Figure 3. Density plot of the magnetic field distribution (15), 
which was chosen for our computations. Parameter values 
here correspond to the most probable case with the central 
field strength B(0, 0) = 50 /iG. For more details see subsection 
lie. Color version of our article is available online. 



model parameter. These alternative less realistic values 
were chosen to be 5(0, 0) = 15 /iG and 5(0, 0) = 300 /iG. 
They can be considered as approximate boundaries of the 
interval where the actual field value lies. The detailed dis- 
cussion of the final results variation due to uncertainties 
in 5(0,0) will be done in the subsection IIIB. 

At this point we have completely specified magnetic 
field distribution. For illustration purposes we showed 
the density plot of our distribution in eq. (15) on fig. 3 
for the case of 5(0, 0) = 50 /iG. 

The magnetic field distribution in eq. (15) is written 
in a cylindrical coordinate system attached to M31. Axes 
x'y^ lie in the M31 disk plane. Another non trivial step 
is to make a transformation between coordinate system 
attached to M31 and "laboratory" system with Oz axis 
pointing from the M31 center to the observer. This is 
necessary for the flux computation, because we are going 
to perform a numerical integration over the emitting vol- 
ume along our line of sight, which is coincident with Oz 
axis and not coincident with Oz'. The two coordinate 
systems are shown on fig. 4. We performed the coor- 
dinate transformation (p^cp^z) < — > {p'^cp'^z') using the 
formulas below, which were derived in [24] (formula (6)) 
for cartesian coordinates (x^y^z) < — > {x' ^y' ^z'): 



x' = —X sinP -\-y cos P, 

y' = —X cos P cos i — y sin P cos i — z sin i, 

z' = —X COS P sin i — y sin P sin i -\- z cos i. 



(16) 



where P and i represent the position angle and axis in- 
clination respectively, and define the orientation of the 
M31 disc plane with respect to the sky plane. Follow- 



ing [6], we assumed: P = 38°, i = 78°. From eq. (16), 
using relations between cylindrical and cartesian coordi- 



nates it's easy to obtain the one-to-one correspondence 

{p,^,z) i — >\p',ip',z')'. 



((— p cos ip sin P -\- p sin cp cos P)^ + (p cos cp cos P cos i -\- p sin cp sin P cos i -\- z sin i) 
—p cos (p cos P sin i — p sin (/:? sin P sin i -\- z cos i. 



2\2 



(17) 



Finally, we substituted eq. (17) into eq. (15) obtaining 
therefore P as a function of {p^ip^z), which allowed to 
compute the final fiux. 



D. Estimated fluxes 

A final expected radio flux due to DM annihilation in 
M31 halo is defined by eq. (1). In order to compute it 
we combined all relevant relations (2)-(4), (6)-(15), (17) 
and substituted them in eq. (1). The integration limits 
for spatial coordinates can be easily figured out from the 
fig. 4: 



^ p ^ {d — z) tana, 

0^(p^27T, 

-rioo ^ ^ ^ rioo, 



(18) 



where a is the angular radius of our ROI. Fig. 5 presents 
the dependence of the total expected radio fiux at z^ = 74 
MHz on the angular radius of ROI for both selected an- 
nihilation channels and three trial WIMP masses: m^ = 
10, 100 and 1000 GeV. Boundaries of every shaded region 
correspond to the two limiting cases of DM density dis- 
tribution in the halo, which were discussed in subsection 
II B: cioo = 12 and cioo = 28. Thus the actual fiux is ex- 
pected to lie somewhere inside shaded regions. The thick 
central lines present the algebraic averages between the 
corresponding boundary curves (fiuxes). The fiux satu- 
rates around a ~ 5' and does not grow significantly at 
larger radii. We also computed the expected signal-to- 
noise ratio (SNR) for the same set of parameters and 
presented it on fig. 6. Details about noise level computa- 
tions can be found in subsection III A. These plots justify 
our choice of a. Basically, in this choice we are governed 
by two counteracting effects. On one hand, fig. 6 clearly 
suggests to use the smallest a possible in order to achieve 
the best SNR. On the other hand, as discussed in subsec- 
tions 2.2-2.3, the accuracy of DM density and magnetic 
field strength estimates decreases when approaching the 
central point of the halo. As a matter of fact, we know 
almost nothing about real distributions of these quanti- 
ties inside the region of size a ^ 0.25'. Hence, our ROI 
should be much larger than ^ 0.25' in radius in order to 
encompass a volume where our overall knowledge of pa- 
rameters distribution is precise enough. Based on these 
considerations, we chose a = 5' as the optimal ROI ra- 
dius. With such choice we lose little in SNR, however. 



I 

we gain in overall reliability of calculations. Also, we can 
not increase a any higher than ^ 5', because it will lead 
to non-desirable capturing of projected point sources. 

As an intermediate step of calculations, we decided to 
study also the contribution of leptons with different en- 
ergies into the total fiux. For this purpose we computed 
the part of the total fiux produced by the leptons with 
initial energies at production below different thresholds 
Em- We spanned the whole possible range of Em val- 
ues from the lepton's rest energy to the WIMP's rest 
energy. The results are presented on fig. 7 for all four 
frequencies used. These plots are generated for the most 
relevant magnetic field model, WIMP mass m^ = 100 
GeV (relevance of this mass scale will be seen in the re- 
sults section) and the Cioo = 12 halo model. Essentially, 
we can see on fig. 7 that the main contribution to the 
total fiux (> 50%) for the all frequencies and annihila- 
tion channels comes from a relatively narrow window of 
initial lepton energies between ^ O.Olm^c^ and O.Itti^c^, 
or 1 GeV and 10 GeV. This energy window of one order of 
magnitude width is considered to be narrow with respect 
to the whole energy range of produced leptons between 
~ 10~^ m^c^ and Im^c^, which spans more than 5 orders 
of magnitude. This implies that the most important frac- 
tion of the lepton population for our considerations has 
an energy above ~ 1 GeV. This fact, in turn, means that 
according to fig. 2 the main energy loss mechanisms for 
such population would be synchrotron and ICS emission. 
Thus, these kinds of losses dominate in comparison with 
the other two - bremsstrahlung and Coloumb losses, and 
play a main role for final results. Meanwhile, it's also 
useful to note that in the case of r~^r~ channel the to- 
tal fiux is constituted by more energetic leptons than the 
fiux for bb channel, which is in complete agreement with 
the difference between these two channels mentioned in 
the Introduction - r~^r~ channel produces harder lepton 
energy spectrum. 

Another important result of our computations is the 
frequency dependence of radio fiux from M31 halo. These 
dependencies are shown on fig. 8 for three different 
WIMP masses, both annihilation channels, both limit- 
ing halo models and the chosen ROI with a = 5^ As a 
general feature, we see, for all cases, a decrease of the ra- 
dio fiux as the frequency increases. The steepness of this 
decrease is, however, more pronounced for lower WIMP 
masses. This behavior suggests to primarily leverage on 
low frequency observations for obtaining DM constraints. 
However, the sensitivity of radio surveys typically in- 




Figure 4. Two coordinate systems {x,y,z) and {x\y'^z'). The plane x'y' is coincident with the M31 disc plane. Integration 
over the emitting volume is performed in {x,y,z). More details are in subsection II C. 
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Figure 5. The expected radio flux dependence on the angular radius of the ROI. The expected flux lies inside shaded regions. 
Upper and lower boundaries of these regions correspond to the different halo concentration values cioo = 28 and cioo = 12 
respectively. Thick central lines present algebraic averages between corresponding limiting cases. Curves are shown for the 
two annihilation channels XX ~^ bb, t~^t~, three WIMP masses m^ = 10,100, 1000 GeV and the most probable magnetic 



field distribution with the central value of 50 /iG. The annihilation cross section is (av) 
frequency is z/ = 74 MHz. For more details see subsection II D. 
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Figure 6. The expected SNR dependence on the angular radius of ROI. The expected SNR lies inside shaded regions. Upper 
and lower boundaries of these regions correspond to different halo concentration values cioo = 28 and cioo = 12 respectively. 



Curves are shown for the two annihilation channels XX 



T^T , three WIMP masses m^ 



10, 100, 1000 GeV and the 



most probable magnetic field distribution with the central value of 50 /xG. The annihilation cross section is {av) =3-10 
cm^/s, the observational frequency is z/ = 74 MHz. For more details see subsection II D. 
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Figure 7. The part of total radio flux produced by leptons with energies at production below different thresholds. These 
illustrations are generated for the following set of parameters: cioo = 12, {cjv) = 3 • 10~^^ cm^/s, B(0, 0) = 50 /iG. WIMP mass 
is chosen to be vn^ = 100 GeV, which reflects the most relevant mass scale with respect to final results, as will be seen below. 



We can note that the majority of the total flux are produced by leptons with initial energies above 
are in subsection II D. 
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Table I. Main parameters of all radio surveys used. 



Survey 


Frequency, 
MHz 


Beam 

diameter 

or FWHM 

arcseconds 


RMS noise 

level inside 

the beam 

a[''\ mJy 


VLSS 


74 


80 


100 


WENSS 


325 


82 


4.0 


NVSS 


1400 


46 


0.5 


GB6 


4850 


240 


4.0 



creases with frequency, partially compensating the effect 
of the lower signal. There is, therefore, no obvious op- 
timal frequency to be used and that would lead to the 
strongest constraints possible. For this reason, we opted 
radio observations at very different frequencies for actual 
constraints derivation. This procedure and its results will 
be described in the next section. 



III. OBTAINING CONSTRAINTS ON DM 
ANNIHILATION 

A. Comparison Avith radio surveys 

In order to make a next step on the way to final con- 
straints, we needed to obtain the upper limits on the 
actually observed radio fluxes from our ROI. Compari- 
son of these upper limits with the expected fluxes com- 
puted above immediately provides the upper limits on the 
annihilation cross section for different values of WIMP 
masses, since the expected fluxes is directly proportional 
to {(Jv). For this purpose we studied many major radio 
surveys conducted in the past and chose four of them 
(corresponding frequencies and detailed descriptions are 
provided in brackets): VLSS (74 MHz, [25]), WENSS 
(325 MHz, [26]), NVSS (1400 MHz, [27]) and GB6 (4850 
MHz, [28]). All other surveys were rejected due to ei- 
ther one of the following two reasons: too low resolution 
of a survey (i.e. the telescope beam size is larger than 
the ROI size) or not covering M31 area of the sky. In 
table I we outlined the main surveys' parameters, which 
were collected from [25]- [28] and then used in calcula- 
tions. The uncertainties of survey frequencies, caused by 
finite bandwidth, do not exceed the level of ~ 4%, and 
do not affect final results significantly. 

Real sky images were viewed and analysed by the Al- 
adin software, which is described in [29]. The images of 
our ROI from all four radio surveys are shown on fig. 



9. Firstly, we notice that our ROI does not capture any 
contaminating point sources. On the last three images 
(325, 1400 and 4850 MHz) some signal from M31 center 
is clearly visible. The 74 MHz image presents essentially 
noise. For the precise flux upper limit derivations we 
needed to obtain the measured signal values q (z = 1 — 4 
denotes different observational frequencies) and the rms 
(root of mean square or root of dispersion) noise values 
(Ji for our ROI for all images (we assume Gaussian distri- 
butions of noise levels inside a radio telescope beam and 
inside the ROI). The measured signal values were easily 
read by the Aladin. As for the rms noise values, we de- 
rived from basic principles that they can be estimated 
as ai = al ^/Wl = cry ^^ where Ni are the numbers of 
beams contained in the ROI, Pi are the angular radii of 
beams (or half of the FWHM for each survey) and cri 
are the rms noises inside a beam for each survey (cited in 
table I). Thus, we have obtained both necessary ingredi- 
ents for a flux upper limit estimation from the ROI - the 
measured signal values and the rms noise values. 

Having all the necessary information, we obtained the 
limiting cross section values for cases of different param- 
eter sets using the following formula for the probability 
density of the noise values inside our ROI: 



Pi{ni 



-Wi{{(Jv))) 



1 



27r<Ji 



exp 



[cj - Wi{{crv))f 
2a? 



(19) 
where n^ - the noise values inside the ROI; q - the mea- 
sured signal values from the ROI and Wi{{av)) - the ex- 
pected signals from WIMP annihilation, which is linearly 
proportional to an unknown annihilation cross section 
{(Jv). Here we do not include any background radiation 
(this question will be further developed in the next sec- 
tion). For our final constraints onto (av) we decided to 
use 99.73% confedence level, which corresponds to 3a 
Gaussian confedence level. 

In order to guide intuition on which experiment is the 
most constraining for a given particle mass, in fig. 10 we 
showed our results on the cross section when each fre- 
quency is considered individually. There we plotted the 
constraints obtained at 99.7% confidence level by four 
different surveys used, for two halo models (described 
in subsection II B), two annihilation channels selected 
XX ~^ ^^7 r~^r~ and the most probable magnetic field dis- 
tribution with 5(0, 0) = 50 /iG. We spanned the range of 
WIMP masses between 6 GeV and 1000 GeV and com- 
puted the limiting (crv) values for a subset of WIMP 
masses inside this range, interpolating results for mass 
values in between. The general key properties of our 
results can be summarized as follows: first of all, differ- 
ent halo models yield very different limiting annihilation 
cross sections (by about one order of magnitude) for all 
WIMP masses. This reflects the high sensitivity of our 
results to the halo model assumed. Another noticeable 
fact is that r~^r~ annihilation channel appears to be less 
constrained than bb channel, in complete agreement with 
the mentioned expectation that harder secondary spec- 
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Figure 8. The frequency dependence of the expected radio flux for three different WIMP masses and two annihilation channels. 
ROI radius a — 5^ shaded regions reflect uncertainty in DM density distribution as on the previous plots, 5(0,0) = 50 /iG, 
{cFv) — 3 • 10~^^ cm'^/s. More details are in subsection II D. 



tra channels produce generally lower final radiation flux. 
It is also interesting to notice that, for different WIMP 
masses and annihilation channels, the most constraining 
survey among all the ones considered varies. Therefore, 
we can not isolate a specific key frequency/experiment 
for each plot, bearing the role of the most constraining 
for all WIMP masses - we need to use all experiments to- 
gether to get the best constraints. And final constraints 
from the joint analysis will be gathered in the next sec- 
tion together with analysis of the role of magnetic field 
uncertainties. 



B. Joint analysis and final constraints 

Constraints shown on fig. 10 present exclusion con- 
tours from four independent observations. As we can 
see, these contours lie rather close to each other on each 
diagram at least for some WIMP masses. Such situa- 
tion suggests to apply a joint likelihood analysis in order 
to infer combined constraints from all four independent 
measurements, which we can expect to be better than the 
constraints taken just considering the exclusion curves in 
fig. 10. 

Signals other than DM are expected at each frequency. 
An appropriate modeling of all other astrophysical sig- 
nals, with their frequency dependence, is necessary in 
the joint analysis of various frequency maps. We de- 
cided to make no assumptions on astrophysical processes 
other than DM annihilation, and simply modeled the to- 
tal contribution of all other background signals as an 



a-priori unconstrained offset uncorrelated between fre- 
quency bands. We then performed a standard Bayesian 
analysis of our data. Following the Bayes theorem (for 
more details see [30]): 



p{Q\c) 



p(^e)p(e) 

p{c) 



(20) 



where the l.h.s. is the desirable probability density dis- 
tribution for unknown set of parameters O in case of the 
observation of data c, p{c\^) = L{S) is the probabil- 
ity density of our observed data c for fixed parameters 
O (or likelihood function), p{Q) is the prior probability 
density of B and p{c) is the marginal likelihood. In our 
case 6 = {{crv) ,5^) is the set of unknown parameters, 
which we are aiming to infer - annihilation cross section 
and background radiation fluxes Si on all four frequen- 
cies involved. c= (q) is the data measured. Then the 
likelihood function L will have the form: 



p{c\e) = L(e) 



JJexp 

i=l 



(Ci 



■w,i{av))r^ 



2a 



(21) 

This likelihood function presents essentially the product 
of the independent noise level Gaussian distributions on 
all frequencies written in terms of the measured and ex- 
pected fluxes (taking into account that Ci = Si-\-Wi-\-ni). 
As for p(B), we made no specific assumptions and chose 
a constant (flat) prior. The function p{c) in eq. (20) does 
not require specification, since it does not depend on 6 
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Figure 9. Radio images of the M31 central region. Our ROI of 5^ radius is marked by the red circle. Corresponding surveys 
are commented on each image. We can find noise only on the 74 MHz and significant signals on the other frequencies. 



and, hence, plays a role of a normalization constant. Af- 
ter we constructed the five- dimensional density distribu- 
tion (20) for all five unknown quantities, we needed to in- 
tegrate eq. (20) over all possible ranges of non-interesting 
parameters (marginalization). In our case we were not 
interested in background values, so we obtained the dis- 
tribution for the desired quantity {av) by marginalizing 



over the background values: 

p{{crv)\c)r^ f L{e)p{e)ds 

' {Ci -Si -Wi{{(Jv))Y 



]^exp 



2(7? 



dsi..ds4,^ 

(22) 



where we took into account the flat prior. Thus formula 
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Figure 10. All the computed constraints: the WIMP annihilation cross section values above the contours are excluded at 99.7% 
confedence level. Constraints from different surveys are marked by different colors and thicknesses, which are explained on one 
of the plots and same for all other plots. Standard thermal reUc value ^gv) — 3 • 10~^^ cm^/s is also shown by the red dashed 
line. Corresponding halo model and annihilation channel are commented on each plot. No emission sources besides DM are 
assumed. Magnetic field has the most probable distribution with B(0,0) = 50 /xG. For more details see subsection III A. 



(22) gives us the final probabiHty density distribution 
for annihilation cross section combined from all four ob- 
servations with possible unknown backgrounds included. 
From this distribution we easily constructed the final 
exclusion contours. We presented them for the case of 
99.73% confidence level on fig. 11-12, where the allowed 
area is below the exclusion contours for the correspond- 
ing models. 



We included in our final results possible uncertainties 
in the magnetic field distribution discussed in subsection 
lie. As we outlined there, this distribution is uncertain 
in two aspects - its vertical scale height z'^ and the cen- 
tral field strength 5(0,0). As for the z'^^ our trial runs 
showed that Zq variation over all possible values is able to 
change the final exclusion values of (av) by no more than 
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{o'v)_2Q% for all possible sets of model parameters. This 
variation is considerably smaller than the one induced 
by, e.g. halo model uncertainties; and would suggest less 
conservative limits. For these reasons, we do not con- 
sider it in our final results. However, we can not treat 
in this way the second important parameter - 5(0,0), 
because practice showed significant variation of the fi- 
nal exclusions due to uncertainties in 5(0,0) discussed 
in subsection II C. According to this discussion, we ex- 
pect 5(0,0) = 50 fiG to be the most probable, but we 
also consider here 5(0, 0) values of ^ 15 and ~ 300 /iG 
as limiting cases. We obtained all our final exclusions for 
these three different magnetic field distributions. 

Our results are presented in fig. 11-12. In these fig- 
ures, the dashed lines represent the specific halo models 
(either cioo = 12 or cioo = 28), while the solid lines corre- 
spond to the algebraic averages between the two limiting 
halo models, for fixed magnetic field values. The two 
figures show limiting curves for different combinations 
of the (three) possible magnetic field models and (two) 
possible halo models. Fig. 11 shows more explicitly the 
dependence on the field strength, while fig. 12 shows 
the dependence on the halo model for the reference field 
strength 5(0,0) = 50 /iG, and allows to compare these 
uncertainties with the ones caused by variation in the 
magnetic field. 

As we can see, annihilation with the standard thermal 
relic cross section (av) = 3- 10~^^ cm^/s is excluded for a 
significant part of WIMP mass range. We listed in table 
II the precise WIMP mass values for all models consid- 
ered, below which (av) is less than 3 -10"^^ cm^/s. These 
masses were obtained as the contact points between the 
corresponding exclusions and the thermal relic threshold 
on fig. 11-12. We can see a significant spread in limiting 
masses around those which correspond to our preferred 
model with 5(0,0) = 50 /iG and the averaged fiux be- 
tween two limiting halo cases. For this most realistic 
model the WIMP masses which do not allow the stan- 
dard relic cross section are the ones smaller than 100 
GeV and 55 GeV for bb and r~^r~ channels respec- 
tively. In general, these limits were derived with rather 
conservative assumptions. Discussion of our results and 
comparison with other studies will be done in the next 
section. 



IV. SUMMARY AND DISCUSSION 

In our work we aimed to obtain robust constraints on 
WIMP-mass/annihilation-cross-section plane by study- 
ing radio observations of M31. Annihilating DM in M31 
halo produces highly energetic secondary leptons, which, 
in turn, emit synchrotron radiation due to Andromeda 
magnetic field. We computed the expected characteris- 
tics of this radiation. We first modeled the DM density 
distribution in M31, using the standard NFW density 
profile. Parameters of this distribution were found in the 
previous studies of the M31 halo. These studies appeared 



to have good agreement between each other in halo mass 
determination (Mioo ^ 10^^ M©), but they differ signif- 
icantly in the obtained halo concentration parameter - 
Cioo ^ 12 vs. Cioo ^ 28. That's why we treated these two 
cases separately and considered them as limiting cases of 
possible halo models for M31. Another important step of 
the calculations was magnetic field distribution modeling. 
For this purpose we used available magnetic field mea- 
surements in M31. Information, which is not available 
for the M31, was gathered by extrapolating MW mag- 
netic field distribution properties and general properties 
of galactic magnetic fields as well. The idea to use MW 
magnetic field measurements was motivated by high sim- 
ilarity between the M31 and the MW, and, also, by the 
much more detailed knowledge we have about the MW. 
We introduced a realistic field distribution, which is char- 
acterized by the 50 /iG value of the central field strength 
and exponential decline both in vertical and radial-in- 
plane directions. Then we computed the secondary par- 
ticle yields from WIMP annihilation using the results of 
DarkSUSY package. We showed that neglecting the spa- 
tial diffusion of the injected leptons is a reasonable ap- 
proximation. We calculated final fiuxes in optically trans- 
parent media approximation. As for the sky area chosen 
for putting the constraints, we chose the circle with the 
angular radius 5' around the M31 center, and motivated 
our choice in subsection II D. After obtaining the radio 
fiuxes expected, we compared them with available radio 
observations of M31. For this purpose we used four radio 
surveys spanning wide range of frequencies: VLSS (74 
MHz), WENSS (325 MHz), NVSS (1400 MHz) and GB6 
(4850 MHz). We analyzed the data from these surveys 
and compared them with the theoretical expectations for 
radio fiuxes. We then derived limits on the cross section 
{(jv) as a function of particle mass (ranging from 6 GeV 
to 1000 GeV). Our analysis is conservative in sense of 
the absence of any specific assumptions about a possible 
radiation other than from DM. The final results are pre- 
sented in the section III on fig. 11-12. We also computed 
our results for the alternative less probable magnetic field 
models with the central field values of 15 and 300 /iG in 
order to study the dependence of our final constraints on 
this model parameter. 

Our main result is the exclusion of the possibility 
of WIMP annihilation with the thermal cross section 
(av) = 3- 10~^^ cm^/s or higher for WIMP masses below 
^ 100 GeV (^ 55 GeV ) for the bb ( r^r') annihila- 
tion channel. However, these estimations are affected by 
significant uncertainties. Taking into account potential 
uncertainties in DM density and magnetic field distri- 
butions, these limits could be as low as 23 GeV and as 
high as 280 GeV for the bb channel. For annihilation into 
r~^r~ these limits respectively are 23 GeV and 130 GeV 
(see table II). 

Our results are comparable with other findings on 
WIMP masses derived through other observations. The 
best up-to-day constraints in indirect DM searches were 
obtained by joint analysis of MW satellites by Fermi 
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Figure 11. Joint constraints including all radio surveys and backgrounds. These plots show all potential uncertainties of both 
distributions of DM density and magnetic field strength. The dashed contours reflect specific halo models discussed in the text 
(Unes represent 99.7% confedence level). The continuous lines present averages between the exclusions for the corresponding 
limiting halo models. Different continuous contours illustrate the variation of exclusions due to magnetic field uncertainties. 
The upper and lower exclusions enclose the shaded regions of all other possible exclusions from the most conservative to the 
most optimistic ones. More details are in subsection IIIB. 
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Figure 12. Joint constraints obtained. These plots incorporate potential uncertainties of both distributions of DM density and 
magnetic field strength. The continuous line presents the most realistic exclusion contour, which corresponds to the result of 
our work on the fiducial case we assumed (average exclusion between those which correspond to the two limiting halo cases and 
the most realistic magnetic field distribution with B(0, 0) = 50 /iG). The shaded area around the continuous contour, which is 
constrained by the inner dashed green lines, corresponds to dark halo model uncertainties only. Larger shaded region, which 
is enclosed into the outer dashed lines, incorporates both uncertainties of dark halo and magnetic field. For more details see 
subsection IIIB. 
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Table II. WIMP masses in GeV, below which annihilation with the canonical cross section is excluded, for the different halo 
and magnetic field models. The most realistic cut off values are emphasized. 



Central field 
B(0,0) 


Cioo = 12, 
XX -^ bb 


Averaged 
fiux, XX -^ bb 


Cioo = 28, 
XX -^ bb 


Cioo = 12, 

XX -^ T+T- 


Averaged fiux, 

XX -^ T + T- 


Cioo = 28, 

XX -^ T^T~ 


15 ^J,G 


23 


45 


160 


23 


35 


88 


50 /LtG 


63 


100 


280 


39 


55 


130 


300 /LtG 


75 


110 


280 


38 


53 


110 



gamma telescope. This work is presented e.g. in [3], 
where the authors obtained the constraints by compari- 
son of Fermi data with the expected gamma emission pro- 
duced directly in WIMP annihilation into the same chan- 
nels. They reported the lowest allowed WIMP masses 40 
GeV and 19 GeV for hh and t^t~ channels respectively. 
When considering 95% systematic errors on the DM dis- 
tribution within the dwarves, their derived WIMP mass 
lower limit varies between 19-240 GeV {hh) and 13-80 
GeV {t^t~). The constraints derived in this work, sug- 
gesting masses above 100 GeV and 55 GeV in order to 
obtain a relic cross section, are within the range of un- 
certainties of current Fermi results. This suggests the 
relevance of multi-frequency approach in indirect detec- 
tion: all wavelength ranges may make valuable contri- 
butions into an overall progress. As for other studies 
dedicated to indirect searches on radio frequencies, we 
can compare our results with the similar studies for MW 
and M33 - [8] and [18] respectively. Our constraints are 
significantly stronger than the last cited, since the ex- 
clusion contours obtained by the authors of [8], [18] lie 
above the thermal relic value {crv) = 3 • 10~^^ cm^/s 
on the WIMP-mass/annihilation-cross-section plane for 
any possible WIMP masses. 

We show here, for the first time, that it is possible 
to constrain the WIMP parameter space by radio 
observations of M31. Specifically, while making con- 
servative assumptions, we were able to exclude small 



WIMP masses for conventional, thermally produced 
DM. However, the range of higher masses of hundreds 
GeV is still absolutely allowed and unexplored. This 
paper presents only the first relatively simple step in a 
comprehensive analysis of M31, and does not yet include 
all potential uncertainties and relevant effects. This 
work may be expanded in two main directions. Firstly, 
astrophysical uncertainties can be treated with a more 
general approach: including, for example, spatial diffu- 
sion of annihilation products and possible variations in 
magnetic field modeling. In addition, while encouraged 
by current finding, the potential to study this object 
with other current and upcoming radio telescope may 
be explored. 
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Appendix A: Absorption analysis 

Here we discuss all potentially relevant mechanisms of 
absorption of radio emission during its propagation from 
the M31 center to the observer. We focused only on a 
potential absorption inside M31 and did not consider pos- 
sible absorption in MW, because the Andromeda galaxy 
lies far enough from the Galactic plane, and we do not 
expect any significant absorption in this direction on the 
sky. In general, inside the Andromeda there are two main 
agents of interstellar medium relevant for our purposes: 
a dust and an ionized gas. Among them a dust can be 
excluded as an absorber rather easily: typically a dust 
grain size does not exceed 10~^ cm (see e.g. [31]). And 
the wavelengths of radio emission involved are 10 — 1000 
cm, which is many orders of magnitude larger than grain 
sizes. In such case, when a wavelength is much greater 
than an obstacle size, it is well known that radiation does 
not interact with obstacles. Thus, we can conclude that a 
dust does not affect radiation in the considered frequency 
range. 

A situation with the absorption by interstellar plasma 
is more tricky. Here we can distinguish several possible 
ways of absorption. First of all, let's check the Langmuir 
frequency of plasma involved. As well known, a radiation 



can not propagate through plasma, if its frequency is less 
than the Langmuir frequency of propagation medium. 
The Langmuir frequency is defined as 



^pi 




(Al) 



where n is the concentration of plasma electrons, which 
needs to be estimated. For this purpose we used the 
results of [10], where the gas density distribution was ob- 
tained for the central region of MW. Taking into account 
that MW and M31 are very similar galaxies, we can ex- 
trapolate results for MW on M31 staying at necessary 
level of accuracy. From [10] we can easily see, that the 
plasma concentration in a galactic center can not exceed 
n < 10 cm~^ in the worst case, which yields Vpi < 30 kHz. 
Thus, the Langmuir frequency is much smaller than the 
observational frequencies of MHz-GHz, which means that 
the interstellar plasma is transparent with respect to this 
absorption mechanism. 

Another potentially relevant mechanism is a syn- 
chrotron self-absorption. Synchrotron emission, gener- 
ated by relativistic leptons from DM annihilation, can 
be absorbed by neighbour emitting leptons. In order to 
estimate the level of this absorption quantitatively we 
computed the corresponding optical depth along our line 
of sight, which goes through the M31 center: 



OLaaQji ^ 



(A2) 



los 



where ags is the synchrotron self- absorption coefficient. 
We used the derived expression for ass from [32] (formula 
(6.50) there): 



87ri/2 



/ dEPe{E,v,r) 



where Pe{E^ v^ f) is the synchrotron power of one lepton 
defined by eq. (3)-(4), ^^{E^f) is the stationary energy 
spectrum of emitting leptons of one kind derived by eq. 
(6). Then we substituted all relevant values of the pa- 
rameters involved. Particularly, for WIMP mass and fre- 
quency of observations we used as an example m^ = 100 
GeV and u = 74 MHz. Such combination of all parame- 
ters yields Tss ~ 10~^, which means that we can neglect 
by synchrotron self-absorption completely and reliably. 
And this conclusion is valid for all other frequencies used 
in our work, because ass decreases with frequency, and 
74 MHz is the lowest frequency of observations. In gen- 
eral, this conclusion agrees with the results of [5], where 
the authors solved similar DM annihilation problem with 
application to MW. And they also showed an irrelevance 
of synchrotron self- absorption on all frequencies. 
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And the last potential absorption mechanism, which 
has to be checked, is a free-free absorption by plasma. 
Here we also estimated the corresponding optical depth: 



estimated as 



Tff = J affdl, 



(A4) 



los 



where the absorption coefficient aff was also taken from 
[32] (formula 5.18b there): 



aff = 3.7 • lO^T-^^^Z^n^u-^ ( 1 



exp 






(A5) 

where n and T are plasma concentration and tempera- 
ture respectively. As for the spatial distribution n{f)^ 
it does not appear to be obtained in the literature for 
M31. That's why we decided to use such distribution 
for MW again as an approximation. Taking all relevant 
information from [10] and assuming gff ^ 1 and Z = 1 
(hydrogen plasma) we estimated Tff < 0.01 in the worst 
case scenario. Thus, we can see that the free- free ab- 
sorption is not relevant for our work as well in the first 
approximation. 

Summarizing, in this section we have conducted the 
detailed analysis of all potentially relevant absorption 
mechanisms of radio emission generated by DM annihila- 
tion. None of these mechanisms achieves significant level. 
This conclusion is in agreement with similar studies [5], 
[8] for the MW. Thus we can ignore any absorption in 
our analysis without significant loss of accuracy. 



Appendix B: Spatial diffusion of annihilation 
products 



Tloss{E) 



E 
E 



(B2) 



_ h{E,ry 
where h{E, f) = hies + hsync + hrem + hcoi (see eq. (8)- 
(11)). After all necessary substitutions we obtained Id 
dependence on the lepton energy E^ which is shown on 
fig. 13. For the magnetic field B and the concentra- 
tion n we used the expected values for the M31 center 
~ 50 jiG and ~ 0.1 cm~^ respectively. These parame- 
ter values were justified in the section II. As we can see 
on fig. 13, the lepton diffusion path over relevant range 
of energies does not exceed ~ 500 pc. Our emitting re- 
gion, which we capture by our ROI with angular radius 
a ^ 5^ would have the form of the cylinder with the 
radius pmax ^ Oid ^ 1100 pc. Thus, the smallest size 
of the emitting region is about 2 times larger than the 
diffusion length of leptons in the M31 center. It means, 
in turn, that the leptons do not have enough time to mi- 
grate significantly and escape from the emitting volume, 
before they cool down and discontinue radiation. Also 
taking into account the fact that the majority of total 
radiation flux due to DM annihilation is formed in the 
very central region, we can conclude that inclusion of the 
spatial diffusion in our calculations should not affect our 
results significantly. And neglecting by the diffusion is an 
acceptable approximation in our computation procedure. 
Such conclusion is in general agreement with the work [8] , 
where the similar procedure of constraints derivation was 
conducted for the MW. However, since the spatial scale 
of the emitting volume and the lepton diffusion path do 
not differ drastically, we allow an opportunity to include 
the diffusion in our calculations in a future work in order 
to improve accuracy of our results. 



Here we study the role of the spatial diffusion of anni- 
hilation products in our problem. In order to understand 
the importance of the spatial diffusion for final results, we 
should compare the characteristic distance, which anni- 
hilation products travel while they are emitting relevant 
radiation, with the characteristic size of the emitting re- 
gion. According to e.g. [8], the diffusion path traveled by 
leptons can be calculated as Id = \/D{E)rioss{E)^ where 
D{E) is the diffusion coefficient for leptons participating 
in eq. (5); rioss{E) is the cooling time of leptons, dur- 
ing which they are emitting the expected radiation and 
losing their energy until escaping the relevant energy do- 
main. We took the diffusion coefficient D{E) from [8], 
where it is provided for the MW, and expected to work 
roughly for the Andromeda galaxy as well: 



I>(^)=I)o(|) , 



(Bl) 



with 1^0 = 10^^ cm^/s, Eq = 3 GeV and a Kolmogorov 
spectrum 6 = 1/3. The cooling time for leptons can be 
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Figure 13. The dependence of the lepton diffusion length Id 
on the energy E. The magnetic field used is B = 50 /xG, 
concentration is n = 0.1 cm""^. 



